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Abstract 

We derive a free energy for weakly anisotropic ferromagnets which is 
valid in the whole range of temperature and interpolates between the mi- 
cromagnetic energy at zero temperature and the Landau free energy near 
the Curie point T c . This free energy takes into account the change of the 
magnetization length due to thermal effects, in particular, in the inho- 
mogeneous states. As an illustration, we study the thermal effect on the 
Stoner- Wohlfarth curve and hysteresis loop of a ferromagnetic nanoparticle 
assuming that it is in a single-domain state. Within this model, the saddle 
point of the particle's free energy, as well as the metastability boundary, 
are due to the change in the magnetization length sufficiently close to T c , as 
opposed to the usual homogeneous rotation process at lower temperatures. 



PACS numbers: 75.10.-b, 75.50.Tt 
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1 Introduction 

The macroscopic energy of weakly anisotropic ferromagnets which first appeared 
in the seminal paper by Landau and Lifshitz JT|, has been an instrument for 
innumerable investigations of domain walls and other inhomogeneous states of 
magnetic systems. Later an approach based on this macroscopic energy has been 
called "micromagnetics" ||. Strictly speaking, micromagnetics is an essentially 
zero-temperature theory for classical magnets, as it considers the magnetization 
as a vector of fixed length. Under these conditions, it can be easily obtained as 
a continuous limit of the classical Hamiltonian on a lattice. Practically, micro- 
magnetics has been applied to nonzero temperatures as well, with temperature- 
dependent equilibrium magnetization and anisotropy constants. 

On the other hand, close to the Curie temperature T c a magnetic free energy 
of the Landau type can be considered. This free energy allows the magnetization 
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to change both in direction and length, and it is, in fact, a continuous limit 
of the free energy following from the mean- field approximation (MFA). Using 
this approach, Bulaevskii and Ginzburg || predicted a phase transition between 
the Bloch walls and the Ising-like walls at temperatures slightly below T c . Much 
later this phase transition was observed experimentally [f|, || using the theoretical 
results for the mobility of domain walls in this regime || |7|. 

Although both of these approaches have been formulated at the same mean- 
field level, they have been considered as unrelated for a long time. On the other 
hand, the MFA itself is an all-temperature approximation, and thus one can 
question its macroscopic limit in the whole temperature range. The corresponding 
macroscopic free energy should interpolate between the Landau-Lifshitz energy, 
or micromagnetics at T = and the Landau free energy in the vicinity of T c . The 
Euler equation for the magnetization which minimizes this generalized free energy 
appears in Refs. and ||. Whereas the condition for the Landau theory is 
M -C M s , where M s is the saturation magnetization at T = 0, the new equations 
only require \M — M e \ <C M s , where M e is the equilibrium magnetization at zero 
field. For weakly anisotropic magnets (the anisotropy energy is much less than 
the homogeneous exchange energy, and this is satisfied by most compounds) the 
magnetization magnitude M is either small or only slightly deviates from M e , 
thus the condition above is satisfied in the whole temperature range. 

Derivation of the generalized macroscopic free energy from the MFA is much 
subtler than that of the Euler equations. This free energy appears in Ref. || 
where its form was guessed. In this paper, we will present this derivation and il- 
lustrate the resulting magnetic free energy in the case of a single-domain magnetic 
particle with a uniaxial anisotropy. For the sake of transparency, we will con- 
sider the exchange anisotropy rather than the single-site anisotropy; this does not 
change the qualitative results. It will be shown that the free-energy landscape for 
the magnetic particle has different forms at low temperatures and near T c . At low 
temperatures, the saddle point between the "up" and "down" minima is located 
near the sphere |M| = M e , the value of M being somewhat reduced in compari- 
son with M e due to thermal effects at T > 0. Near T c , the value of M strongly 
changes by going from one minimum to another through the saddle point; for zero 
fields the saddle point becomes M = 0. These effects also modify the metasta- 
bility boundary of the magnetic particle (the famous Stoner-Wohlfarth curve 
which was experimentally observed in Ref. fH|), and its hysteresis curves. 

The main body of this paper is organized as follows. In Sec. |2| we give the 
derivation of the magnetic free energy at all temperatures. In Sec. |3| the free- 
energy landscape of a single-domain magnetic particle is analyzed. In Sees. f| 
and [5] we study the Stoner-Wohlfarth curves and hysteresis loops at different 
temperatures. In Sec. |B| we discuss the possibility of observing the new thermal 
effects for single-domain magnetic particles. 
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2 Free energy of a magnetic particle 



Let us start with the biaxial ferromagnetic model described by the classical 
anisotropic Hamiltonian of the type 

H A^O-H ^""^ — ^ ' Jiji^SziSzj ~\~ Vy^yiSyj Vx^xi^xj) i (1) 

i ij 

where //q is the magnetic moment of the atom, i,j are lattice sites, Sj is the 
normalised vector, |s$| = 1, and the dimensionless anisotropy factors satisfy r\ x < 

T}y < 1. 

To study the macroscopic properties of this system at nonzero temperatures, 
it is convenient to use a macroscopic free energy. In the literature one can find 
two types of macroscopic free energies for magnets. One of them is the so-called 
micromagnetic free energy which is valid at zero temperature, and the other one is 
Landau's free energy which is applied near the Curie temperature T c . Since both 
free energies are based on the mean- field approximation (MFA), it is possible to 
derive a simple form of the MFA free energy for weakly anisotropic ferromagnets 
which is valid in the whole temperature range and bridges these two well-known 
forms. 

The free energy F = —T In Z of a spin system described by the Hamiltonian 
in Eq. ([!]) can be calculated in the mean-field approximation by considering each 
spin on a site % as an isolated spin in the effective field containing contributions 
determined by the mean values of the neighboring ones. Namely, 

n n MFA = n 00 -J2 H f FA ^ ( 2 ) 

i 

where 

^00 = 9 X! "A? ( a zi<^zj + VxVxiVxj + VyVyi^yj) > (3) 

(Ti = (sj) is the spin polarization, and the molecular field H^ FA is given by 

Hf FA = /i H + J ij {°zje z + r] x a xj e x + r]y^yj e y) ■ ( 4 ) 

3 

Then the solution of the one-spin problem in Eq. (|2]) leads to 

F = H 00 - NT ln(47r) - T ]T A(&) 

i 

A(f)Eh(^), (5) 

where N is the total number of spins, = and ^ = /5H^ FA . The MFA free 
energy determined by Eq. (Q) and Eq. (|5|) can be minimized with respect to the 
spin averages erj to find the equilibrium solution in the general case where the 
anisotropy 1 — r) x>y is not necessarily small. The minimum condition for the free 
energy, dF/d(Ti = 0, leads to the Curie- Weiss equation 

<T i = B{£ i )^, (6) 

Si 
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where B(£) = coth(£) — l/£ is the Langevin function. 

For small anisotropy, 1 — T] X}y <C 1, one can go over to the continuum limit 
and write for the short-range interaction Jy 

JijVj = Jocri + J aA(Ti, (7) 

3 

where A is the Laplace operator acting on the components of <x(r), J is the 
zero Fourier component (the zeroth moment), and JqO is the second moment of 
the exchange interaction Jy. For the simple cubic lattice with nearest neighbor 
interactions a = a^/z, z = 6, and is the lattice spacing. Going from summation 
to integration in Eq. (El), one obtains 



F NT 

Jq J( 



T . . If fl 2 1 . , 1 
- ln(47r) + - dr\ -a 2 + -(a, h eS - h) - — 

VqJ [Z Z (JJq 



A(0 , (8) 



where vq is the unit-cell volume, 

£ = (3J {ct + h eff ) 

h cff = h + aAcr - (1 - rj x )a x e x - (1 - r] y )a y e y , (9) 

and h = /x H/J . We will consider the case of small fields, h -C 1. Since in this 
case in Eq. (0) |h cff | <C |<x| in the whole range below T c (near T c the value of a 
is small but the susceptibility is large), the last term of Eq. (g) can be expanded 
to first order in h eff using 

£ = 6 + <J£, £o = AV, <^/3J — ^ (10) 

(7 

and the first two terms of the expansion 

A(0 = A(Co) + BfcOtf + ^'(&)(^) 2 , (11) 
where = dB(g)/d£. Hence, 

= — / rfr i o a ~~ «r A ^° CTh ~ o CT ' hcff ~ h r 

J V (2 (JJq a \ a 2 J J 

NT 

ln(4vr). (12) 



J, 



o 



Near T c = T C MFA = Jo/3 the order parameter er becomes small, and using A(£) = 
£ 2 /6 — ^ 4 /180 and B(£ ) = £ /3 = /3 c J cr/3 = o one straightforwardly arrives at 
the Landau free energy 



J = ^/*{^(v-) 2 --h + 5 (1 -"^ + 5 (1 -">»-f' T2 + l c 

NT 

- ^-ln(4vr), (13) 

•A) 

where e = (T C MFA - T)/T C MFA and (W) 2 = (VaJ 2 + (V^) 2 + (Va,) 2 . Note 
that Eq. ( |T3"D formally yields unlimitedly increasing values of o at equilibrium, 
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as a function of the field h. To comply with the condition a <C 1 in Landau's 
formalism, h should be kept small, as was required above. One could also work 
out the h 2 corrections to Eq. ([13]). 

Now we consider the temperature region where the influence of anisotropy and 
field on the magnitude of the spin polarization a can be studied perturbatively. 
Concerning the influence of anisotropy, the applicability criterion can be obtained 
from the requirement that the local shift of T c for spins forced perpendicularly to 
the easy axis should be smaller than the distance from T c , i.e., AT C /T C ~ 1 — i] <C 
e. The macroscopic free energy in the perturbative region can then be combined 
with the Landau free energy in their common applicability range 1 — i] <C e <C 1. 
Thus, in the perturbative region we expand the first two terms of expression Eq. 
(O) up to the second order in 5a = a — o e using 



6) = & + Ze = PJ Ve, 6^/3J 5a (14) 

and the formula analogous to Eq. (|TT|). This leads to 



^- = ^ + -fdr (1(1 - B'(3J )(a - a e f -^A*h - ( ^1 _ 1 ) {(7 . h „. _ h) y { 15) 
J J v J 12 a \ a 



where B' = B'(£ e ), and 

F e NT 



ln(47r) + iV 
■Jo Jo 



(16) 



is the equilibrium free energy in the absence of magnetic field and the quantity 
a e is the spin polarisation at equilibrium satisfying the homogeneous Curie- Weiss 
equation 

*e = B{Q. (17) 

The minimum condition for Eq. (P~5[) , 5F/5er = 0, after an accurate calcula- 
tion taking into account the dependence of -B(£o) on c an d neglecting the terms 
quadratic in h e g, results in an equation of the form 

j-( g - g .)--b«+ |gx[ ';_ xM1 =o, us) 

XII o- a 2 x\\ 

where h e ff is given by Eq. (|9]) and the dimensionless longitudinal susceptibility x\\ 
is given in Eq. (|T9|) below. The solution of Eq. (|l^) satisfies cr||h e fr, and the term 
with the double vector product plays no role. Considering the response to small 
fields h = h z e z and h = h x y e x ^ y in Eq. (|T8| ) in a homogeneous situation (in the 
transverse case h c s = and a = cr e ), one can identify the reduced susceptibilities 
for the spin polarization as 

da z B'/3J 

X\\ ~ 



dh z l-B'flJo 
1 _ 1 

1 - V* 1 - Vy 

Our expression for the free energy, Eq. (|i"5|) , is still cumbersome, but can be 
simplified if we make the observation that in the perturbative region the deviation 
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8a = a — a e is proportional to h e $, and the terms of the type 5a • h cS and {8a) 2 in 
Eq. ([H]) are thus quadratic in h e s. Such terms are nonessential in the calculation 
of F itself, they are only needed for the proper writing of the equilibrium equation 
Eq. (|i8|). Now we can replace Eq. ( |T5D by a simplified form 

T = f + vj dl {^ (Vct)2 - CTh + 1 (1 - "> 2 - ">» + ' (20) 



This form coincides with Eq. (|15|), if we set a = a e , i.e., it yields the same value 
of F in the leading first order in h e g. On the other hand, Eq. (^) leads to the 
same equilibrium equation, Eq. ([18]), without the nonessential last term. Now, 
at the last step of the derivation, one can combine Eq. (|20|) with the Landau free 
energy in Eq. (Jf3|), leading to 

y = T + 1 M U«(V-) 2 - -h + ^-a\ +-L^ + -4-(a 2 - a 2 ) 2 ) . (21) 
J J V [2 2xx 2x?; 8<7|X|| J 

Indeed, in the Landau region, t < 1, from Eq. (|i"7|) and Eq. (]TQ|) it follows 
a e 2 = (5/3)e and X|| = (2c)" 1 , and Eq. (glj) simplifies to Eq. In terms of the 

magnetization M defined by 

M(r) = »o<T{r)/v , *{r) = (s r ), (22) 

and other dimensional quantities, the free energy Eq. (^l|) takes on the form (5|. § 

F = F - + ft { i (VM)2 " MH + + 4 M » + s^ (M2 - M - 2)2 } m 

with 

9d = —r, Xa = —rXa, W D = , (24) 

where is the so-called dipolar wave number, Wd is the characteristic energy 
of the dipole-dipole interaction, and Xa = dM a /dH a with a = x,y,z are the 
susceptibilities [cf. Eq. flTtp]. The applicability of Eq. (|3|) requires that the 
deviation M — M e from the equilibrium magnetization M e in the absence of 
anisotropy and field, is small in comparison with the saturation value M s = fio/vo- 
This is satisfied in the whole range of temperature if the anisotropy and field are 
small, i.e., 1 — r) X) y <C 1 and [i^ll <C Jo. On the other hand, the free energy in 
Eq. (^) can be transformed into the "micromagnetic" form by introducing the 
magnetization direction vector v = M/M. One can then write 

' < = V5, K *.v = t¥-, (25) 



^Ax,y ^AX,y 

where K XjV are the anisotropy constants. In particular, for the uniaxial model 
one can rewrite 

7^-(M 2 + M 2 ) = -Kv 2 z +K, (26) 
2x± 

where K is the uniaxial anisotropy constant. Note that at nonzero temperatures 
M and thus the anisotropy constants can be spatially inhomogeneous. In this 
case Eq. (E3T) is more useful than its micromagnetic form. 
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It should be stressed that the traditional way of writing the magnetic free 
energy in terms of the magnetization is, at least from the theoretical point of 
view, somewhat artificial. All terms in Eq. (|23"D apart from the Zeeman term, are 
non- magnetic and in fact independent of the atomic magnetic moment //q. The 
latter cancels out of the resulting formulae, as soon as they are reexpressed in 
terms of the original Hamiltonian (]l|). On the other hand, Eq. fl2"3f) is convenient 
if the parameters are taken from experiments. 

Examination of the formalism above shows that it can be easily generalized to 
quantum systems, leading to the same form as in Eq. fl2~3"|). In the derivation, the 
classical Langevin function B(x) is replaced by the quantum Brillouin function 
B s (x). 



3 The free-energy landscape for a uniaxial mag- 
netic particle 

Henceforth, we will consider the uniaxial anisotropy, r\ x — rj y = rj±. For single- 
domain magnetic particles in a homogeneous state, the gradient terms in the free 
energy can be dropped and the free energy of Eq. (|23| ) can be presented in the 
form 

F = F e + (VM 2 / X ±)f 

f = -n-h+ 1 -(n 2 x + n 2 y ) + ±-(n 2 -l) 2 J (27) 

where V is the particle's volume and / the reduced free energy written in terms 
of the reduced variables 

n = M/M e , h = U X ±/M e , a = 2 X \\/x±- (28) 

One can see that the parameter a here controls the rigidity of the magnetization 
vector; it goes to zero in the zero-temperature limit (the fixed magnetization 
length) and diverges at T c as a = (1 — rj)/e within the MFA. As the MFA is not 
quantitatively accurate, it is better to consider the susceptibilities and hence a 
as taken from experiments. Although this procedure is not rigorously justified, it 
can improve the results. Note that the reduced free energy / in Eq. (p7|) is only 
defined for T < T c since the reduced magnetization n is normalized by M e . 

For fields h inside the Stoner-Wohlfarth astroid, which will be generalized 
here to nonzero temperatures, / has two minima separated by a barrier. Owing 
to the axial symmetry, one can set n y = for the investigation of the free energy 
landscape. The minima, saddle points, and the maximum can be found from the 
equations df /dn x = df /dn z = 0, or, explicitly 

n z {n 2 — 1) = ah z 

n x (n 2 — 1 + a) = ah x . (29) 
These equations can be rewritten as 

1 - n 2 = -ah z /n z = a(l - h x /n x ), (30) 
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whereupon the important relation follows 



hx h z h x Tl z /o-i\ 

= 1 or n x = . (31) 

n x n z h z + n z 

Using the latter, one can solve for n x and obtain a 5th-order equation for n z 

h\n\ = (h z + n z ) 2 (ah z + n z -n 3 z ), (32) 

from which the parameters of the potential landscape such as energy minima, 
saddle point, and energy barriers can be found. 

In zero field, the characteristic points of the energy landscape can be simply 
found from Eqs. (p9|). One of these points is n x = n z = 0, which is a local 
maximum for a < 1 and a saddle point for a > 1. The minima are given by 
n x = 0, n z = ±1. The saddle points correspond to n z = 0, while from the second 
of Eqs. fl29|) one finds 



J ±Vl - a, a < 1 , . 

n * = { 0, a > 1. ( 33 ^ 

In fact, due to the axial symmetry, for a < 1 one has a saddle circle n x +n 2 = 1 — a 
rather than two saddle points. The free-energy barrier following from this solution 
is given by 

Af = f - f - / ( 2 ~ a )/ 4 ' a ^ 1 (Wi 
J - Jsad Jmin \ l/(4a), a>l. 1 J 

The free-energy landscape in zero field is shown in Fig. [[]. At nonzero temper- 
atures a > 0, the magnitude of the magnetization at the saddle is smaller than 
unity since it is directed perpendicularly to the easy axis, and for this orientation 
the "equilibrium" magnetization is smaller than in the direction along the z axis. 
For a > 1, the two saddle points, or rather the saddle circle, degenerate into a 
single saddle point at n x = n z = 0, and the local maximum there disappears. 
That is, for the magnetization to overcome the barrier, it is easier to change its 
magnitude than its direction. This is a phenomenon of the same kind as the phase 
transition in ferromagnets between the Ising-like domain walls in the vicinity of 
T c (the magnetization changes its magnitude and is everywhere directed along 
the z axis) and the Bloch walls at lower temperatures 0, f3|. 



4 The Stoner-Wohlfarth curve 

The Stoner-Wohlfarth curve separates the regions where there are two minima 
and one minimum of the free energy. On this curve the metastable minimum 
merges with the saddle point and loses its local stability. The corresponding 
condition is 

d 2 f/dn 2 x x d 2 f/dn 2 z - (d 2 f /dn x dn z ) 2 = 0, (35) 

or, explicitly 

(1 - n 2 ) (l - 3n 2 - a) + 2an 2 z = 0. (36) 



S 




Figure 1: The free energy of a ferromagnetic particle with uniaxial anisotropy 
[/ in Eq. (p7|) in zero field] for a = 2x\\/x±- — 0-5 (upper plot) and a = 2 (lower 
plot) corresponding to lower and higher temperatures, respectively. 

Using Eq. (|3"0"|), one can transform the equation above to the quartic equation for 

n z 

h z [(2 + a)n z + 3ah z ] + 2^ = 0. (37) 

Before considering the general case, let us analyze the limiting cases a <C 1 and 
a > 1. 

At low temperatures, i.e., a < 1, the magnetization only slightly deviates 
from its equilibrium value, and from Eq. (136|) , to first order in a, one obtains 

n 2 = 1 - an\ or n\ + (1 + a)n\ = 1. (38) 

From Eq. ( |3~7"|) and the analogous equation for n x and using Eq. (|30|) , one can 
derive the field dependence of n z and n x on the Stoner-Wohlfarth curve 

n z = -hl/ 3 [l-(a/2)(h^-l/3)} 

n x ^h^[l + {a/2){hT-l)\. (39) 

Inserting these results in Eq. ([$]), one arrives at the equation for the Stoner- 
Wohlfarth astroid 

h 2 J 3 + [(l + a/2)h z f 3 ^l, a<l, (40) 
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Figure 2: Dependence h z (a) at h x = on the Stoner-Wohlfarth curve. 

where a is given in Eq. (pS| ). One can see that, in comparison with the standard 
zero-temperature Stoner-Wohlfarth astroid, i.e. at a = 0, h z is rescaled. The 
critical field in the z direction decreases because of the field dependence of the 
magnetization magnitude at nonzero temperatures. 



In the case a ^> 1, i.e. near T c , Eq. (j36|) relating n x and n z in the equation 
for Stoner-Wohlfarth curve simplifies to 

n 2 x + 3n 2 z = l. (41) 

Using this equation together with Eqs. (|30|), one obtains 

n z = -(ah z /2) l '\ n x S h x . (42) 

Making use of this result in Eq. ([|I]), one obtains another limiting case of the 
Stoner-Wohlfarth curve 

hi + 3(ah z /2) 2/3 = 1, a > 1. (43) 

In this case, the critical field (i.e., the field on the Stoner-Wohlfarth curve) in 
the z direction is strongly reduced, and there is no singularity in the dependence 
h cz (h x ) at h x = 0. 

The qualitatively different character of the Stoner-Wohlfarth curves in these 
two cases is due to the different mechanisms pertaining to the loss of the local 
stability for the field applied along the z axis. For h x = the mixed derivative 
(<9 2 / / dn x dn z ) in Eq. ( p5| ) vanishes, and Eq. fl35|) factorizes. Explicitly, we have 

(a-l+n 2 z )(-l + 3n 2 z ) = 0. (44) 

Vanishing of the first factor in this equation corresponds to the loss of stability 
with respect to the rotation of the magnetization, d 2 f /dn 2 = 0. Vanishing of 
the second factor, d 2 f/dn 2 = 0, implies the loss of stability with respect to the 
change of the magnetization length. Using n x = 0, with the help of the first 
equality in Eqs. fl30D one obtains in the two cases 




_ . _= VI -a, a < 2/3 , . 

hcz ~^ »• -2/(3 3 / 2 a), a>2/3. 1 
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Figure 3: The Stoner-Wohlfarth curves at different temperatures, a = 2x\\/X-i 
(T = 0), 0.3, 0.5, 2/3, 0.835, 1, 2. 



Note that the transition between the two regimes occurs here at a different value 
of a than in Eq. (^). The dependence h z (a) at h x = is shown in Fig. |2|. 

In the general case, it is easier to find the Stoner-Wohlfarth curve numerically 
from Eqs. ( p9|) and fl3"6"|). The results in the whole range of a are shown in Fig. [3[ 



5 Hysteresis loops 



In this section we use our model to study the hysteresis loop, i.e., the dependence 
n(h) at the stable or metastable free-energy minimum. At first we consider the 
case h x = in which the problem can be solved analytically. Setting n x = 



in the first of Eqs. (|29|) one obtains the cubic equation n z 



ah 2 



0, the 



solution of which, for the positive branch of the hysteresis curve, reads 




.1/3 



+ 



ah, 



1/3 



\h z \ < h z 
\h z \ > h. 



(46) 



where h z \\ is given by Eq. fl45|) and 

<p = &TCCOs(h z /h z n), D = (a/2) 2 (h 2 z - h 2 zll ). 



(47) 



The negative branch of the hysteresis curve can be obtained by the reflection 
h z =^ —h z and n z =^ — n z . Eq. (|46|) is written in the form which is explicitly real. 
In fact, both forms hold in the whole range h z > —h z \\, and there is no change of 
behavior at h z = h z n, as expected on physical grounds. The characteristic values 
of the magnetization on the positive branch are n z (— h z \\) = 1/V3 and n z (h z \\) = 
2 / y^. The solution above was obtained by setting n x = and thus ignoring 
the transverse instability which can occur before the longitudinal instability at 
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Figure 4: Thermal effect on the static hysteresis loops [a(0) = 0, a(T c ) = oo]. 
Upper plot: tj) — (h||e 2 ), lower plot: ip = tt/4- 

/i 2 = — for the positive branch. This competition of instabilities has been 
studied in the previous section. It was found that for a < 2/3 the transverse 
instability occurs at the fields \h z \ = \h z \\\ = a/1 — a < \h z ±\. Thus in this case 
the branches of the hysteresis curves should be cut at ±/i*n; at these fields the 
system jumps to the other branch. 

The analytical results obtained above for h x = and the numerical ones for 
the case of h directed at the angle ip = tt/A to the easy axis are shown in Fig. 
|j. For h x = 0, the derivative dn z /dh z diverges at the metastability boundary for 
a > 2/3 (longitudinal instability). 

6 The homogeneity criterion, thermal activa- 
tion, spin-waves 

In the preceding sections we illustrated how the general magnetic free energy 
of Eq. (|23| ) works for the simplest model of a uniaxial magnetic particle in a 
single-domain state at elevated temperatures. Here we discuss the possibility 
of observing the new types of behavior found above. This requires satisfying 
two rather restrictive conditions. First, the particle spends in the metastable 
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minimum a time that is long enough for measuring only if the barrier height 
energy is much larger than thermal energy At any T > 0, the particle will escape 
from the metastable state via thermal activation with a rate exponentially small 
at low temperatures. For this reason, strictly speaking, the static hysteresis at 
nonzero temperatures does not exist and dynamic measurements are needed. At 
elevated temperatures, the required frequency of these measurements can become 
too large. Second, the free energy barrier in the single-domain state should be 
smaller than the barrier energy related with the formation of a domain wall which 
would travel through the particle and switch the magnetization from one state 
to the other. The two criteria can be combined as follows 

T < AF SD < AF DW . (48) 

The first criterion here requires that the particle's volume is high enough whereas 
the second criterion requires that it does not exceed some maximal value. Let 
us consider, for instance, the zero field case, in which Fsd is the single-domain 
free-energy barrier of Eqs. ([34]) and (p7|) . At not too high temperatures a<l, 
one has Af ~ 1, and the domain walls are usual Bloch walls with the energy per 
unit area 

w = 2M 2 J(q 2 d S), S=^xZ/q d , (49) 

where 5 is the domain-wall width. For a spherical particle of radius R, the saddle 
point of the energy corresponds to the domain wall through the center of the 
particle. Using the definitions of parameters introduced in Sec. El, one obtains 
the formula 

AF SD Rq d r R 



, . v 1 - V—, ( 50 ) 

where clq is the lattice spacing. It is seen that the single-domain behavior of 
particles with R^> requires small values of the anisotropy 1 —rj. Working out 
the first inequality in Eq. (0), one can rewrite these equations in the form 

( 9 \ 1/3 R^ 1 

<— ^^^=, (51) 



\(1 - rfial) ^ a ~ y/1 - 7) 

where 6 = 7yT c MFA and a e is the spin polarization at equilibrium. Clearly, for 
^ < 1 these conditions can be satisfied. Observing the new effects exhibited by 
the hysteresis suggested above requires a = 2x\\/x±. = 2x\\/x± ~ 1> which for 
small anisotropy, i.e. 1 — rj -C 1, requires approaching T c . Indeed, near T c the 
first of Eqs. ( |I9"D yields x\\ — (2e) _1 , thus a ~ 1 implying e ~ 1 — 77. In this region 
in Eq. (|5T|) one has 9 ~ 1 and a 1 = (5/3)e. Then the existence of the interval for 
R/a Q in Eq. (|5"lf) requires (1 — rj)/t e which is impossible since in this region 
(1 — 77) /e ~ 1 and e < 1. A similar analysis shows that also in the region where 
a > 1 Eqs. flilf) cannot be satisfied. 

Thus we are led to the conclusion that the qualitatively different types of 
behavior of single-domain magnetic particles for a > 1 cannot be observed with 
the standard techniques. If the particle's size is small enough, the single-domain 
criterion is satisfied, but increasing temperature to a ~ 1 causes strong thermal 
fluctuations. The potential landscape shown in Fig. [TJ is still valid but for studying 
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the dynamics of the magnetic particle in this range we need a special kind of 
Fokker-Planck equation for non-rigid magnetic moments. Such an equation has 
not been considered yet. 

On the other hand, for large particle sizes the energy barriers are high and 
thus the process of thermal activation is suppressed, which is favorable for the 
observation of hysteresis loops. In this case, however, the barrier states are 
those with a domain wall across the particle. Analyzing these states goes beyond 
the scope of this article. We only mention that for a > 1/2 the structure of 
domain walls in a ferromagnet is completely different from that of a Bloch wall: 
The transverse magnetization component in the wall is zero everywhere whereas 
the longitudinal component changes its magnitude and goes through zero in the 
center of the wall 0. For the observation of thermal effects in the hysteresis via 
inhomogeneous states, higher values of the anisotropy 1 — n are needed. In this 
case, thermal effects manifest themselves starting from low temperatures. An 
example of a strongly anisotropic material is Co, for which one obtains 1 — rj ~ 
0.02. This value is still much smaller than one, so that the validity condition for 
the magnetic free energy of Eq. ( P3|) is satisfied. 

Returning to the results of this paper, we can say that one can only observe 
corrections to the well-known results, such as that of Eq. (fEJ), i n the range 
a < 1, i.e., not close to T c . On the other hand, one should not forget that the 
mean-field approximation used in this paper, while leading to qualitatively correct 
predictions, fails to account for some subtler effects that can be responsible for 
important modifications of the results. One of these effects is the influence of spin 
waves on the longitudinal susceptibility which enters the definition a = 2x\\/x±- 
Whereas within the MFA x\\ rapidly decreases with decreasing temperature below 
T c and is independent of the anisotropy, it is infinite in the whole range below 
T c for isotropic ferromagnets because of spin waves. The square-root singularity 
in the dependence M(H) at zero field has been experimentally observed and 
reported on in Ref. ||11|| . In uniaxial ferromagnets, x\\ becomes large for small 
anisotropies. This means that if the mean-field expression for a is replaced by 
its value taken from experiments, which is much larger due to spin-wave effects, 
the thermal effects discussed in this paper will considerably increase in intensity. 
Such a redefinition of a is, of course, not rigorously justified, although it captures 
the essential physics. A more involved approach taking into account spin-wave 
effects for an exactly solvable model confirms the concomitant increase of thermal 
effects on the variation of the magnetization length, as was explicitly shown for 



domain walls in Ref. [TI| 



Finally, we would like to mention that quite recently, experimental results 
have been obtained by Wernsdorfer et al. [jn| on 3 nm cobalt nanoparticles which 
clearly show the disappearance of the singularity near H x = at a temperature 
circa 8 K (the blocking temperature being 14 K). The height of the experimental 
astroid decreases nearly as its width with increasing temperature, but it does not 
become flat as predicted by our calculations, which is not surprising considering 
the fact that T -C T c , but the disappearance of the singularity is definitive. 
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7 Conclusion 



In this paper we have derived a macroscopic free energy for weakly anisotropic 
ferromagnets which is based on the mean-field approximation and is valid in the 
whole range of temperature interpolating between the micromagnetic energy at 
T = and the Landau free energy near T c . As an illustration, we have con- 
sidered single-domain magnetic particles with uniaxial anisotropy and we have 
shown that thermal effects qualitatively change the free-energy landscape at suf- 
ficiently high temperatures, so that the passage from one free-energy minimum 
to the other is realized by the uniform change of the magnetization length rather 
than the uniform rotation. This also qualitatively changes the character of the 
Stoner-Wohlfarth curve and hysteresis loops. The latter effects cannot be ob- 
served with standard methods, however, because keeping the height of the free- 
energy barrier much larger than thermal energy requires so large particle sizes 
that the single-domain criterion is no longer satisfied. For the uniform states, the 
theory is valid at low temperatures, but then the thermal effects considered in 
the paper are small corrections to the zero-temperature results. Large thermal 
effects on hysteresis should be searched for at temperatures close to T c in parti- 
cles of larger sizes, where the saddle point of the free energy is an inhomogeneous 
state. Investigation of the corresponding more complicated processes is beyond 
the scope of this paper. 
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